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The cluster of galaxies A754 undergoes a merger of several large structural units. X-ray observations show a nonequi- 
librium state of the central part of the cluster, in which a cloud of cold plasma 500 kpc in size was identified amid the 
hotter cluster gas. The X-ray image of A754 exhibits a brightness discontinuity, which can be interpreted as a shock wave 
in front of a moving cloud of dense gas. The shock parameters are determined from the jump in intergalactic gas density 
using the ROSAT image. The estimated Mach number is Mj = 1.711q 24 at a 68% confidence level. 



INTRODUCTION 

Clusters of galaxies are the largest gravitationally bound 
objects in the Universe. Mergers of clusters result in the re- 
lease of potential energy into the intergalactic medium at a 
level of 10^^ erg, mainly in the form of heating by shock 
waves. Therefore, observations of shock waves in the in- 
tergalactic medium are particularly interesting because they 
provide information for studying the physical processes that 
accompany the cluster formation. 

In this paper we study A754, a rich cluster of galaxies 
at z = 0.0541 in the stage of violent formation. ASCA 
observations of this cluster clearly indicate that the central 
part of the cluster is in a nonequilibrium state and hot and 
cold regions are identifiable in it (Henriksen and Markevitch 
1996). The pattern probably represents the motion of a cold 
dense cloud through a hotter ambient gas. Structures of this 
kind were observed by the Chandra satellite in several clus- 
ters (Markevitch et al. 2000; Vikhlinin et al. 2001; Mazzota 
et al. 2001). The merging parts of the clusters are generally 
expected to move at a supersonic speed, which leads to the 
formation of shock waves. 

The ROSAT X-ray image of A754 exhibits a surface 
brightness discontinuity, which can be interpreted as a shock 
wave in front of a moving cloud of dense gas. The amplitude 
of the surface brightness discontinuity makes it possible to 
determine the jump in gas density at the shock front and, 
hence, the Mach number M — v/vs, where Vj is the shock 
velocity and Vj is the speed of sound in the medium. Unfor- 
tunately, the ROSAT energy range (0.5-2.5 keV) does not al- 
low us to determine the gas temperature at the shock bound- 
ary; i.e., we cannot independently verify the interpretation 
of the brightness discontinuity as a shock wave using the 
temperature jump. 

We quote the physical length scales assuming Hq = 
50 km s^' Mpc^'. The derive shock velocity and Mach 
number are independent of Hq. 



DATA ANALYSIS 

The ROSAT image of the cluster of galaxies A754 is 
shown in Figure 1. The region of cold gas appears as a 
bright cloud at the image center The shape and location 
of a brightness discontinuity to the left from the center are 
consistent with those of a shock in front of a cloud of cold 
gas moving at a supersonic speed. 

Quantitative information about the gas density distribu- 
tion can be obtained by analyzing the surface brightness pro- 
file. The shape of the front of the putative shock is well de- 
scribed by a circle with a radius of 775 kpc. We use the cen- 
ter of this circle as the reference point for radial distances. 
The surface brightness profile was extracted in a ±15° sec- 
tor with respect to the front symmetry axis^ The profile was 
measured in concentric rings of equal logarithmic width; the 
ratio of the inner and outer ring radii was 1.01. We verified 
that varying the reference point changed the measured jump 
in gas density only slightly. 

Measuring the gas density distribution at large radii re- 
quires careful background subtraction and corrections for 
the telescope vignetting. We used S. Snowden's software 
(Snowden et al. 1994) for the data reduction. The soft- 
ware allows one to eliminate the periods of an anomalously 
high background of particles and scattered solar X-ray emis- 
sion. In addition, the exposure maps are computed and all of 
the known background components are fitted and subtracted 
from the image. The final result is the image in a given en- 
ergy range that contains only the cluster radiation, point-like 
X-ray sources, and the cosmic X-ray background. An opti- 
mal ratio of the cluster surface brightness to the background 
level is achieved by using data above 0.7 keV. 

When measuring the X-ray brightness profile, we elimi- 
nated all of the detectable point-like sources and all the ex- 
tended sources unrelated to the main cluster emission. The 
cluster makes a significant contribution to the total bright- 
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The sector opening angle was determined by the size of the brightness 
discontinuity region; the brightness profile for the outer part of the shock 
was obtained in a wider sector (±20°) to improve statistics. 



425 



426 



KRIVONOS ET AL. 



1-1' 



J 

■I -I 



1 filpB t 

1 ' " " t. ^ 



t f-l ' IL L. 



tr cold cloud 




Fig. 1. The ROSAT X-ray image of A754 in the 0.7-2.0 keV energy band. The brightness discontinuity to the left from the central body most likely 
corresponds to the shock. 



ness even at large distances from the center (~ 2 Mpc). 
Therefore, it is hard to identify a region that could be used to 
directly determine the background level. Instead, we derived 
the background level by fitting the brightness profile at large 
distances with the model consisting of a j3 -model (Cavaliere 
and Fusco-Femiano 1976) representing the cluster emission 
and a constant component representing background: 



/ ,N-3j3 + l/2 

/(r)=/o(l + (r/r,)2j +/, 



(1) 



where I{r) is the surface brightness, 4 is the background 
intensity, and is the core radius. 

FITTING THE BRIGHTNESS PROFILE 

The main parameter of the shock motion is the Mach 
number, M = v/v^.The Mach number can be expressed in 
a standard way (Landau and Lifshitz 1988) in terms of the 
jump in any of the gas parameters (density, pressure, tem- 
perature) at the shock front. The most easily measurable 
quantity is the density jump, because the gas density is re- 
lated to the emissivity by a simple relation: e p^. Conse- 



quently, the problem reduces to deprojecting the brightness 
discontinuity. 

The brightness profile can be best deprojected by fitting 
the data by some analytic dependence. Since the outer part 
of the cluster is most likely to have not yet been perturbed 
by the shock, we can assume that the gas density profile is 
described by the standard J3 -model. Inside the shock front, 
we will be concerned with a narrow range of radii in which 
the density profile can be assumed, to sufficient accuracy, 
to be a power-law function of radius. Thus, we have the 
following model of the plasma emissivity: 



sir): 



-3J3 



r>R 



(2) 



where R is the front position, £2 is the post-shock emis- 
sivity, and is the central emissivity for the j3 model. The 
core radius rc was fixed at a typical value of 250 kpc (Jones 
andForman 1984). 

The model emissivity profile was numerically integrated 
along the line of sight and folded with the angular resolu- 
tion of the telescope (^ 25"). Note that the normalization of 
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the P model in equation (|2} can be expressed in terms of the 
pre-shock emissivity, ei . To determine the parameters of the 
above model, we minimized the value for four parame- 
ters (£2,e2/ei,j3,a). Figure |2] shows the derived model of 
the gas emissivity distribution (b) and the surface brightness 
profile (a). 

RESULTS 

The emissivity jump obtained by fitting the brightness 
profile is trivially transformed to the density jump (£2/61 ): 



pi 



j+0.44 
*-0.32' 



(3) 



where the uncertainties are at a 68% confidence level for 
one interesting parameter From the P2/P1 ratio, we find the 
Mach number of the shock with respect to the ambient gas: 



P2 
Pi 



2+(7-1)m2 



(4) 



L71+°'^^ata68% 1 = L71+0.45 



where 7 = 5/3 is the adiabatic index for monatomic gas 
From (0), we obtain Mi 
-0.24 confidence level 68% (see Fig|3} 

The derived value of Mi corresponds to the gas tempera- 
ture jump T2/T1 ~ L72^Q 2g. In principle, measuring such 
a temperature discontinuity would be an independent test of 
the validity of the interpretation of the observed structure 
as a shock wave. Unfortunately, the available Chandra and 
XMM observations of A754 do not allow us to determine the 
gas temperature outside the shock because of the low signal- 
to-noise ratio. Inside the shock, the temperature is measured 
reliably, T2 = 10 keV (Markevitch et al. 2003) from Chandra 
data. 

A useful quantity is the shock velocity relative to the gas 
on the inner side of the discontinuity: 
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M2 = 



f 2 + (7-l)M- 



V27M2-(7-1; 



1/2 



0.66 ±0.07 



(5) 



The measured value of T2 allows us to determine the speed 
of sound and, hence, the absolute shock velocity: V2 = 
1070 ±115 km s-i, vi = 2100;^^{J km i"'. 

The derived Mach number of a likely shock in A754 is 
close to the theoretically expected values for merging clus- 
ters: M = 2 — 3 (Sarazin 2002). Comparison with other 
clusters is hard because observations of merger shocks in 
clusters are is still rare. In two other best-studied cases, 
(lE0657-56-Markevitchet al. 2002, Cyg A-Markevitch et 
al. 1999), high Mach numbers are also inferred. However, 
in A3667 one observes a merger with M ^ 1 (Vikhlinin et 
al. 2001). Clearly, detailed studies of a larger sample of 
merging clusters would be valuable. 
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Fig. 2. (a) The surface brightness profile in the shock region. The distances are measured from the center of curvature of the front along the shock 
propagation; the solid line represents the fitting curve obtained by integrating the plasma emissivity model (b) along the line of sight. 
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Fig. 3. The 68% and 90% confidence regions for two parameters of our model of the intergalactic gas radiation: the gas density ratio at the shock boundary 
and the parameter j3 in the outer part of the cluster [see Eq. |2}]. (a) The cross marks the point that con'esponds to the best- fit parameters; (b) versus 
Mach number. 
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